Spectra of heavy-light mesons at finite temperature 



D.U. Matrasulov* Kh.T. Butanov, Kh.Yu. Rakhimov 

Heat Physics Department of the Uzbek Academy of Sciences, 
0> ; 28 Katartal St., 100135 Tashkent, Uzbekistan 

O ■ and 
^ : F.C. Khanna 

(D ■ Physics Department University of Alberta 

m 



Edmonton Alberta, T6G 2J1 Canada 



^\ Abstract 

' Spectra of heavy-light meson are studied using potential model and thermofield 

i-^i dynamics prescription. The mass spectra of different heavy-light mesons are calcu- 

lated at different temperatures and compared with those at T = 0. It is found that 
^ the binding mass of heavy- light meson decreases as temperature increases. 

0\ . 

1 Introduction 

(N 

^ ; Recent experiments on heavy-light mesons at BABAR [1], CLEO [2] and BELLE [3] 

Q\ \ collaborations stimulated extensive theoretical studies. Within QCD, relativistic potential 

O ■ models, chiral perturbation theory and non-relativistic potential models that are analogs 

of the hydrogen-like atoms. All numerical work is based on lattice QCD gives probably 
^ \ the best results [4]- [10]. 

An effective method to explore hadron properties at finite temperatures is the 
lattice QCD approach that leads to results of higher precision for spectra and transition 
rates. However it requires powerful computer resources. It is therefore interesting to 
explore methods that are simpler and yield interesting results for the hadron spectra 
at finite temperature. Recently such a method has been used to calculate quarkonium 
spectra at finite temperature in potential models [11]. In this paper we extend this 
method to heavy-light mesons with its dynamics described by the Dirac equation. The 
finite temperature effects are included by using thermofield dynamics (TFD), a real time 
finite temperature field theory. We apply the prescription to the Dirac Hamiltonian with 
Coulomb and a linear potential. The physical content of TFD will be described below. It 
is important to note that TFD has been used to study the Casimir effect and quantum 
chaos in Yang-Mills-Higgs system at finite temperature [12]. 

With the results emerging from RHIC (Relativistic Heavy Ion Collider) and the 
anticipated results from LHC (Large Hadron Collider) properties at finite temperature 
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will be needed. In the process of converting baryonic system to quark-gluon plasma and 
subsequently the hadronization of this plasma would have mesons and baryons at finite 
temperature. This would require how these hadrons behave at high temperatures. The 
question of reactions, change of coupling constants and decay rates finite temperatures 
have been considered already [13] . The present calculation provides additional information 
as to how the mass changes with temperature. 

The plan of the paper is to consider briefly the description of the heavy-light 
mesons using the Dirac equation at zero temperature. In section 3, we give a short review 
of TFD in order to make clear the finite temperature formulation in general, stating 
clearly the underlying details of the theory. In section 4, TFD is apphed to calculate the 
spectra of heavy-light mesons at finite temperature. Some concluding remarks are given 
in the final section. 



2 Dirac equation for heavy-light meson 

Starting point for our study of finite-temperature mass spectra of heavy-light meson is 
the Dirac equation written as [14, 15] 

i- = H^ (1) 

where H = dp + P (nig + Ar) — — + Vq. We note that the confining part of the potential is 
included into this equation as a Lorentz scalar (not as fourth component of the 4- vector). 
Introducing the following time scaling r = r~H, we can rewrite this equation as 

i— = rH^. (2) 

OT 

Using standard substitution \E'(r, 6, r) = e~*^'^^(r, 6) and separating angular variables we 
get the time-independent radial equation: 

= rHr'if, (3) 

where 

' [ £ + 7 -m,-Xr-Vo-f) 

being the radial Dirac Hamiltonian. For our purpose we need to represent this Hamilto- 
nian in terms of annihilation and creation operators which are defined as 

1 a! 1 . 1 d 1 

It is easy to see that these operators obey the following commutation relation 

[a, a^] = 1. 

The radial Dirac Hamiltonian (4) in terms of these operators can be written as 
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where 

.l = ^^^(a + at) + ^(a + atf-Z, 

B = -^(a + a^)(a - a^) + 
C = ^(a + a^)(a-a^) + K, 

and 

L> = ^-^(a + aT) - -(a + aT) - Z. 

The eigenvalues of the Hamiltonian written in terms of annihilation and creation operators 
can be calculated by diagonalizing it in a harmonic oscillator basis. The wave function is 
expanded in terms of nonrelativistic harmonic oscillator wave functions as 

Then the energy eigenvalues of the heavy-light meson are to be calculated by diagonal- 
ization of the matrix 

E={^' \ rHr I 

This method will be combined with the thermofield dynamics to calculate spec- 
tra of heavy-light meson at finite temperature. The next section presents briefiy basic 
thermofield dynamics. 



3 Thermofield Dynamics 

Thermofield dynamics is a real-time operator formalism of quantum field theory at finite 
temperature. It has been recognized that for a theory at finite temperature, the Hilbert 
space has to be doubled [16, 17]. The usual Hilbert space may be defined by creation, 
a^{k) and annihilation operators, a{k). Then the second Hilbert space is defined by tilde 
operators, d^{k) for creation and d{k) for annihilation operators. These operators, for a 
boson system, satisfy 

[a{k),a\k')] = Sk,k', 
[d{k),d\k')] = Sk,k', 
[a{k),a'{k')]^0 

and all other commutators are zero. For fermion systems anti-commutator are to be used. 

At finite temperature a BogolyTibov transformation is used to mix these two sets 
of operators to obtain the finite temperature creation and annihilation operators. Any 
operator A{x) is then written at finite temperature as 

A{x,e) ^u{e) A{x) u-\e) 

where U (6) is the Bogolyubov transformation, i.e. 

U{0) = e-'^(^(^» 
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where G — i9{P){a^d^ — aa) with (3 — Here T is the temperature and /cb is the 

Boltzmann constant. 

This imphes that 

a{k, 9) = U{e)a{k)U-\e) = u{e)a{k) - ^{e)a^{k) 

a\k, 6) = U{e)a\k)U-\e) = u{e)a\k) - ^{e)a{k) 

a{k, 9) = U{e)a{k)U-\e) = u{e)a{k) - ^{e)a\k) 

a\k, e) = U{9)a\k)U-\9) = u{e)a\k) - ^{e)a{k) 

where u{e) = cosh(^), ^{9) = sinh(^) such that v?{e) - ^"^{e) = 1. The vacuum state at 
finite temperature is defined similars as 

\o{e)) = u{e)\o,o). 

Then the temperature dependent creation and annihilation operators acting on the vac- 
uum state, I 0(6*)) give 

a{k,9) I 0(^)) = h{k,9) I 0(^)) = 0. 
Thus there is a procedure to work with the two set of operators 

{cAi + Aj) = c*Ai + Aj 
(A) = A 

and [ij, Aj] = 0. It is to be noted that i9^{e) = sinh^ 9 ^ [ef^ - 1]"^ where l3 = 

Prom a group theoretical approach, it has been established that for any kinemat- 
ical symmetry there are two sets of operators, physical observable O and generators of 
symmetry, = — 0. This establishes a connection to an algebraic approach based on 
kinematic symmetry and the doubling does not happen due to an arbitrary ansatz [18]. 

Such a procedure has been used to calculate scattering cross section and decay rates 
at finite temperature. The use of the Bogolyubov transformations that mixes the operators 
from the two Hilbert spaces is a reminder of superconductivity where this implied a 
condensation of quanta in the vacuum state. 

Thus TFD is a powerful tool for exploring quantum dynamics of a system at 
finite temperature provided its Hamiltonian can be represented in terms of annihilation 
and creation operators. It has found many applications in condensed matter physics 
[19, 20, 21], especially in superconductivity theory and related topics. Recently TFD 
has been applied to explore quantum chaos in the Yang-Mills-Higgs system [22] and for 
the calculation of the spectra of a strongly interacting bound system [11]. In this work 
we apply the TFD prescription to the Jaynes-Cummings model. It should be noted that 
TFD has been apphed earlier to the Jaynes-Cummings model [23] where the thermal noise 
effects in quantum optics are studied. In the next section we apply the TFD prescription 
for the calculation of the energy spectra of heavy-light mesons. 
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4 Spectra of heavy-light meson at finite temperature 



Applying TFD prescription to the Hamiltonian of the heavy-hght meson we can write 
it in temperature-dependent form. Then finite-temperature energy eigenvalues of heavy- 
light meson can be calculated by diagonalization of the following temperature-dependent 
matrix: 

{13) bn'AP) 



where 



hn',n{(3) = -\Mf3) + ('^ + ^) ^n',n. 
Cn',n{(3) = ^AsiP) + {k-^) Sn',n, 



and 



Ai{(3) = (x/n 5n',n-i + Vn + 1 Sn',n+i) (cosh 6 + sinh 9) , 

A2{p) = (^^n{n-l) Sn',n-2 + ^(n + l)(n + 2) Sn',n+2 + (1 + 2n)Sn',n^ (cosh^ 9 + sinh^ 9) 

+2 cosh 6^ sinh 6l((n -Fl) 



MP) = [\/n{n - 1) 6n',n-2 - V (n + l)(n + 2) 6n',n+2j (cosh^ 9 + sinh' 9). 

Using this method we can calculate finite-temperature energy and mass spectra of 
heavy-light mesons. 



5 Results 

In table 1 mass spectra of cs-meson is presented for different values of principle and 
orbital quantum numbers at different temperatures. The parameters are chosen as follows: 
rris = 0.5 GeV, rric = 1.486 GeV, A = 0.8 GeV^, Vo = -0.6 GeV, a, = 0.32 {Z = |q;J. It 
is clear from this table that increasing the temperature leads to decreasing of the mass. 
Such effect has been found earlier in the case of quarkonium spectra at finite temperature 
[11, 24]. Table 2 represents the spectrum of 6s-system for the following set of parameters: 
m, = 0.5 GeV, = 4.88 GeV, A = 0.8 GeV^, Vb = -0.6 GeV, a, = 0.22. Table 3 
represents the spectrum of fcc-system for the following parameters: rric — 1.486 GeV, 
rub = 4.88 GeV, A = 0.2 GeV^ Vq = -0.3 GeV, a, = 0.22. 

Again, decreasing of the mass by increasing the temperature is clear from these 
tables. This effect is shown in Fig. 1 which is the pictoral presentation of some states for 
cs-system at finite temperature for the values of temperature 0, 0.1, 0.15, 0.25 (in GeV). 
Results for (bs) and (be) have a similar behaviour of changing with temperature. 
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6 Conclusions 



In this work we have studied spectra of heavy-hght mesons at finte temperature. By com- 
bining quark potential model and real-time finite-temperature field theory, the thermofield 
dynamics Dirac Hamiltonian for heavy-light meson is presented in the temperature- 
dependent form. Diagonalizing this Hamiltonian in the non-relativistic harmonic oscilla- 
tor basis finite-temperature energy and mass spectra of heavy-light mesons with different 
heavy and light quark contents are obtained. It is found that the masses of heavy-light 
mesons at finite-temperature are less than those at T = 0. 

It is important to remark that TFD has proved to be a convenient approach for 
such a study that calculates the effect of temperature on the spectra of heavy-light quark 
mesons [18]. Furthermore the results will provide a useful input to studies of quark-gluon 
plasma at RHIC and LHC. 
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Table 1: The mass spectra (in GeV) of (cs)-meson at finite temperature (in GeV). The 
parameters are chosen as follows: mj = 0.5 GeV, rric — 1.486 GeV, A = 0.8 GeV^, 
Vo = -0.6 GeV, as = 0.32. 



n{l)/T 





0.1 


0.15 


0.25 


1(1= 


:0) 


1.7575 


1.7199 


1.6809 


1.6623 


2 




1.9462 


1.8730 


1.8336 


1.8157 


3 




2.1005 


2.0794 


2.0323 


2.0151 


4 




2.3061 


2.2585 


2.2108 


2.1968 


5 




2.4911 


2.4693 


2.4139 


2.3996 


6 




2.7056 


2.6681 


2.6124 


2.6018 


7 




2.9113 


2.8835 


2.8196 


2.8095 




= 1) 


2.1138 


2.0231 


1.9846 


1.9645 


2 




2.2549 


2.1606 


2.1207 


2.1001 


3 




2.3950 


2.3295 


2.2823 


2.2626 


4 




2.5635 


2.4896 


2.4429 


2.4251 


5 




2.7252 


2.6737 


2.6196 


2.6035 


6 




2.9122 


2.8554 


2.8011 


2.7869 


7 




3.0982 


3.0517 


2.9900 


2.9786 


1(1= 


=2) 


2.5330 


2.4093 


2.3732 


2.3531 


2 




2.6414 


2.5287 


2.4886 


2.4663 


3 




2.7646 


2.6673 


2.6213 


2.6003 


4 




2.9053 


2.8103 


2.7629 


2.7416 


5 




3.0498 


2.9685 


2.9160 


2.8976 


6 




3.2101 


3.1318 


3.0777 


3.0597 


7 




3.3757 


3.3059 


3.2470 


3.2332 
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Table 2: The mass spectra (in GeV) of (6s)-meson at finite temperature (in GeV). The 
parameters are chosen as follows: mj = 0.5 GeV, nib — 4.88 GeV, A = 0.8 GeV^, 
Vo = -0.6 GeV, as = 0.22. 



n{l)/T 









0.1 


0.15 


0.25 


1(1= 


:0) 


5. 


,1481 


4, 


.9984 


4. 


,9756 


4, 


.9317 


2 




5, 


,3432 


5, 


.1995 


5. 


,1741 


5, 


,1230 


3 




5. 


,4986 


5, 


.3576 


5. 


,3361 


5, 


.2764 


4 




5. 


,7082 


5, 


.5695 


5. 


,5445 


5, 


.4757 


5 




5. 


,8936 


5, 


.7548 


5. 


,7343 


5, 


.6574 


6 




6. 


,1098 


5, 


.9715 


5. 


,9476 


5, 


.8603 


7 




6. 


,3158 


6, 


.1767 


6, 


,1585 


6, 


.0625 




= 1) 


5, 


,5057 


5, 


.3519 


5. 


,3269 


5, 


.2898 


2 




5. 


,6503 


5, 


.4992 


5. 


,4712 


5, 


.4251 


3 




5. 


,7907 


5, 


.6414 


5. 


,6165 


5, 


.5607 


4 




5. 


,9627 


5, 


.8159 


5. 


,7877 


5, 


.7233 


5 




6. 


,1253 


5, 


.9817 


5. 


,9591 


5, 


.8858 


6 




6. 


,3147 


6 


.1718 


6. 


,1461 


6 


.0641 


7 




6. 


,5013 


6, 


.3601 


6. 


,3400 


6, 


.2476 


1(1= 


=2) 


5. 


,9259 


5, 


.7708 


5. 


,7454 


5, 


.7147 


2 




6. 


,0364 


5, 


.8820 


5. 


,8535 


5, 


.8138 


3 




6. 


,1600 


6. 


.0062 


5. 


,9780 


5, 


.9270 


4 




6. 


,3028 


6 


.1501 


6. 


,1200 


6, 


.0609 


5 




6. 


,4480 


6, 


.2986 


6. 


,2718 


6, 


.2023 


6 




6. 


,6106 


6, 


.4623 


6. 


,4348 


6, 


.3583 


7 




6. 


,7771 


6, 


.6320 


6. 


,6082 


6, 


.5203 
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Table 3: The mass spectra (in GeV) of (6c)-meson at finite temperature (in GeV). The 
parameters are chosen as follows: = 1.486 GeV, = 4.88 GeV, A = 0.2 GeV^, 
Vq = -0.3 GeV, as = 0.22. 



n{l)/T 





0.1 


0.15 


0.25 


1(1= 


:0) 


6.4132 


6.0423 


5.9721 


5.7302 


2 




6.5940 


6.3623 


6.1911 


5.9703 


3 




7.0008 


6.7803 


6.3437 


6.3006 


4 




7.4043 


7.1883 


6.7436 


6.4567 


5 




7.8131 


7.5976 


7.1403 


6.7371 


6 




7.9541 


7.7572 


7.6095 


7.1157 


7 




8.3235 


8.1021 


7.9986 


7.5723 




= 1) 


7.1444 


7.0931 


6.4040 


5.9279 


2 




7.3652 


7.2417 


7.1465 


6.4397 


3 




7.7834 


7.3741 


7.3373 


7.1281 


4 




7.9654 


7.6891 


7.6270 


7.3249 


5 




8.2946 


7.9966 


7.6666 


7.6077 


6 




8.6206 


8.3817 


7.9385 


7.8986 


7 




9.0263 


8.7682 


8.3063 


7.9818 


1(1= 


=2) 


7.3533 


7.3150 


7.1480 


7.0523 


2 




8.2777 


8.2737 


8.0734 


8.0376 


3 




9.0333 


8.7053 


8.6196 


8.5776 


4 




9.2232 


8.8580 


8.7972 


8.7778 


5 




9.4450 


9.1152 


9.0458 


9.0084 


6 




9.6959 


9.3808 


9.3020 


9.2559 


7 




10.0100 


9.7292 


9.6397 


9.5752 
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Figure 1: The mass spectra of (cs) -meson at finite temperature. 
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